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Abstract 

We investigate the structure of nematic liquid crystal thin films de- 
scribed by the Landau-de Gennes tensor-valued order parameter with 
Dirichlct boundary conditions of nonzero degree. We prove that as 
the elasticity constant goes to zero a limiting uniaxial texture forms 
with disclination lines corresponding to a finite number of defects, all 
of degree \ or all of degree — \ . We also state a result on the limit- 
ing behavior of minimizers of the Chern-Simons-Higgs model without 
magnetic field that follows from a similar proof. 
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1 Introduction 



We investigate disclination line defects in a nematic liquid crystal by using a 
tensor-valued order parameter description based on the Landau-de Gennes 
theory. The unknown field Q in this theory is 5-valued such that Q = 
Q(x,y), where S is the space of 3 x 3, real symmetric, traceless matrices, 
and (x,y) varies in a bounded domain in M 2 . For simplicity, we assume 
that Q is a simply connected bounded domain with a C 3 boundary in the 
plane, representing the reference configuration of a very thin liquid crystal 
material. 

The Landau-de Gennes model is based on a phenomenological theory in 
which stable states of the liquid material correspond to minimizers (or sta- 
ble states) of an energy formulated in terms of Q on f2. The matrix Q(x) 
models the second moments of the orientations of the rod-like liquid crystal 
molecules near x. Its values describe the average orientation and phase of 
the liquid crystals near x, measured through its eigenvectors and eigenval- 
ues. (See Section 1.1 for more detail on this structure.) As such Q is a 
measure of the microscopic anisotropy of their relative positions. In this pa- 
per, we consider fields Q G W 1,2 (Vi; S) with fixed uniaxial nematic boundary 
conditions of the form Q = Qq on <9f2 (in the sense of trace). We assume 
throughout the paper that (Qo)ij € C 3 (dQ) for all 1 < i,j < 3, and 

(1.1) Qo(x, y) = s(n (x, y) <g> n (ar, y) - - I) for (x,y) G dtt 

where / is the 3x3 identity matrix, s is an arbitrary fixed nonzero real num- 
ber, and no is a fixed vector field defined on d£l satisfying no = (ni,7i2,0), 
| no | = 1, and (1.1) on d£l. Note that Qo is invariant under changes in 
direction: no(x,y) — > — no(x,y) at any point (x,y) in dfl,, which allows 
boundary conditions of degree one-half, or integer multiples of one-half, for 
Qq. Nonzero boundary conditions of this type are observed in thin liq- 
uid crystal materials exhibiting defects along curves, known as "disclination 
lines," whose cross-sections in Q are isolated points. (See Figure 1.) We 
analyze a class of equilibria for the Landau-de Gennes energy 



where e > 0, defined for all Q G W 1,2 (Q,S). Here f e is the elastic energy 
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density in Q given by 



fe(Q) = 
+ 

where each term above is summed over all k from 1 to 3, Qij, a denotes 
j^r~, and (xi,X2,x 3 ) = (x,y,z). The above formula is valid in two or 
three-dimensional reference domains. Since here we are considering a two- 
dimensional reference domain S7, we identify Q(x,y) with Q(x,y, 0) above, 
so that Qik,3 = for all 1 < i,j < 3. We assume throughout the paper that 

(1.2) Li > and L x + L 2 + L 3 > 0. 

The term fa is the bulk energy density given by a real-valued C°° function 
which depends on temperature as well as on Q. We assume that temperature 
is fixed and fb = fb(Q) is a nonnegative C°° function defined on S such that 
fb(Q) = if and only if Q G A s = {Q € 5 : Q = s(m m — |7) for some 
m G S 2 } where s is the fixed nonzero constant in the definition of Qq. From 
our definitions in the next subsection, we shall see that the energy well, A s , 
corresponds to a set of uniaxial states. Liquid crystals satisfy the principle 
of frame indifference and are macroscopically isotropic. As a consequence, 
fb is assumed to be invariant with respect to orthogonal transformations, 
that is, we require 

(1.3) fbiRQR 1 ) = fb(Q) for all R G 0(3) and Q G S. 
Set 

S = {Q eS: Q i3 = Q 3i = fori = 1,2}, 



2 Qij,k Qij,k "I - 2 QijJ Qik,k 



L, 



Qij,k Qik,jj 



A = {Q(x,y) G W^in-So): Q = Q on dn}, 

and 

A = {Q G W 1>2 (fi;S): Q = Q on dn}. 

Our goal in this paper is to investigate minimizers for F e in Aq, and to 
analyze their behavior in the vanishing elastic energy limit, e — > 0. The rel- 
evance for doing this is that due to the symmetries described above, these 
minimizers are critical points (equilibria) for the energy F £ over the larger 
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space A, and thus satisfy the full set of Euler-Lagrange equations with re- 
spect to variations in A. (We prove this in Lemma 2.1.) In addition, each 
Q £ S is described in terms of an orthonormal set of eigenvectors. (See 
(1.8).) For Q £ Sq, we have 

(1.4) Q = sim (g> m + 82m 1 - <g> nr 1 - - (s x + s 2 )I 

for some real numbers s\ and S2, and Q has an orthonormal basis of eigen- 
vectors of the form 

(1.5) {m,m^,e3} where |m| = 1, m = (mi, 772,2, 0), 

and rrr 1 - = (-m2,mi,0), 

with eigenvalues 

(1.6) Ai = -(2si - 8 2 ), A 2 = ^(2s 2 - si), A 3 = + s 2 ). 

(See |MN] .) Thus the minimization problem of F e over „4o models the be- 
havior of a thin liquid crystal material occupying f2 x (—rj, rj) with its top 
and bottom surfaces treated so as to fix e3 as a principal axis (eigenvector 
of Q) of the liquid crystal molecules throughout the body, with the other 
two principal axes (eigenvectors) in M 2 x {0}, and boundary values on its 
side given by Q = Qo(x,y). The above problem includes a classic example 
from the liquid crystal literature, in which 

Oh r 

(1-7) / 6 (Q) = / 6 °(Q) = atr(Q 2 )--tr(Q i ) + -(tr(Q 2 )f+D 

= a(E^)-fEAf) + |EA 2 ) 2 +^. 

i=l i=l i=l 

Indeed, taking b, c > 0, a < and an appropriate choice of D, we have 
ft > and / 6 °(Q) = if and only if Q G A s where s = £(b + Vb 2 -24ac). 
(See PN| .) 

1.1 Definitions and Structural Assumptions 

Our results require some structural assumptions on the bulk energy density 
ft). In this section, we state these assumptions, along with some definitions 
and a change of variables in Aq that will be needed to state our main results. 
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It is well known (see |MN] ) that each Q G S has an orthonormal set of 
eigenvectors and can be written as 

(1.8) Q = sin® n + s 2 k ® k- - (s 1 + s 2 )I 

where n and k are orthogonal unit vectors in IR 3 ; moreover, the eigenvalues 
of Q are given by the formula in (1.6). 

Definition 1 Let Q G S. We say that Q is isotropic if all its eigenvalues 
are equal. (In this case, the structure of Q is that of a "normal" liquid.) 

We say that Q is uniaxial if exactly two of its eigenvalues are equal. (In this 
case, Q has an axis of symmetry and its structure is "rod-like" or "disk- 
like".) 

We say that Q is biaxial if all its eigenvalues are distinct. (In this case, 
there is no axis of complete rotational symmetry for Q and its structure is 
"board-like".) 

By formula (1.6) for the eigenvalues of Q G S, it follows that Q is isotropic 
if and only if s± = s 2 = (and hence all eigenvalues are zero); Q is uniaxial 
if and only if one of the following three conditions hold: s\ = and s 2 ^ 0, 
s 2 = and s\ ^ 0, or s\ = s 2 ^ (and hence all eigenvalues are nonzero 
and exactly two of the eigenvalues are equal). Finally, Q is biaxial for all 
other values of si and s 2 . 

The above definition, when applied to a minimizer Q e (x) of F £ in A or Ao, 
allows one to identify subregions of f2 in which the liquid crystal material is 
in an isotropic, uniaxial, or biaxial phase. Note that A s n5o is a disconnected 
set of uniaxial states in Sq with two connected components: A s n So = A' s U 
{s(e3 ® e3 — gj)} where A' s = {s(m (g)m— |/) : m = (mi, m 2 , 0), |m| = 1}; 
also, the boundary values Qo(x,y) are valued in A' s . 

Definition 2 Let 7: [0, 1] — > dQ. be a C 3 positively oriented parameteriza- 
tion of dO, such that 7 is one-to-one on [0,1). For Qq as assumed above, 
choose a unit vector field no(x) = (hi(x),n 2 (x),0) defined on d£l satisfying 
(1.1) such that rio(7(-)) G C 1 ([0, 1)). We define the degree of Qo on d£l by 

— n ( 7 (t)) 1 • " dt: = deg Q . 
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Since limno(7(i)) = ±no(7(0)) by (1.1) and the continuity of Qo, it follows 

that deg Qo = § for some k G Z. Since we are interested in boundary con- 
ditions that correspond to a liquid crystal with dislination-line type defects, 
we assume that k is nonzero, and thus without loss of generality, we shall 
assume throught the paper that k > 0. As e 4- the effect of the bulk energy 
density fj, becomes more pronounced and minimizers tend to have their val- 
ues located in a neighborhood of A s f]So. Due to the boundary conditions, 
however, this cannot happen throughout S7. We prove that the regions in 
which minimizers, Q e (x), of F e take values outside a neighborhood of A' s 
concentrate and quantize into k small subdomains. For a subsequence as 
Sj — > these subdomains tend to k distinct points {a±, . . . , a^} representing 
the cross sections of the limiting disclination lines. 

In |SS| Schopohl and Sluckin carried out a numerical investigation of equilib- 
ria for F £ in A. Their goal was to give evidence that equilibria are strongly 
biaxial near defects. They pointed out that there is a subclass of equilibria 
which is contained in .Ao, an d they developed simulations for these. This is 
the class of solutions that we are studying here. 

To state our main results, we will need the following linear change of vari- 
ables for the coefficients of each Q G Ao in terms of unique functions 
P = (pi,P2) and r: 

Pi + | P2 
Vi | - Pi o 
-r 

From (1.1) and (1.4) each Q G Ao corresponds to a unique (p, r) G W 1,2 (Q; M 2 ) 
W 1 ' 2 ^) satisfying p\ m = p , r\ m = r where |p | = ^, r = §, and 
deg po = k = 2 deg Qo. This can be seen by writing (since no <8> no = 
(-n ) (-n )) 



n o(l(t)) = ±(cos a(t), sina(t), 0) 

for each t in [0, 1), where (cos a(t), sina(i), 0) = no(j(t)) and a G C 1 ([0, 1)). 
Then using (1.1) and(1.4) we observe that 

Po(7(«)) = |(cos2a(i),sin2a(t)). 

We may then recast our minimum problem by considering the set 

A = {(p,r) G WA 2 (fi;R 2 ) x W 1,2 (fi): p = Po and r = S on 30} . 

3 



(1.9) Q = Q(p,r) 
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The mapping Q = Q(p,r) : Ao — > Ao is one-to-one and onto, and the 
eigenvalues for Q(p,r) are Ai = ^ + |p|, A2 = \ — |p|, A3 = — r. By 
(1.3) fb depends only on the invariants of Q; since trQ = 0, these are det 
Q = (IPl 2 - t> and \Q\ 2 = 2|p| 2 + § r 2 . Thus f b (Q) = g b (\p\ 2 ,r) for 
some function g b . We prove in Section 2 that minimizing F £ (Q) over A$ is 
equivalent to minimizing 

(1.10) G £ (p,r) = [ [<? e (Vp, Vr) + £ ~ 2 g b (\p\ 2 , r)} for (p,r) G A), 
where <7 e (Vp, Vr) is defined by 

(l.H) 9e = (L 1 + ^^)|V P | 2 + ( ^ + ^±^)|Vr| 2 
+ ^ 2 1 L3 ^ (p^rj; - pi^ry + r x p2j/ + ^P2x) 

+ |L 2 + L 3 | (pi x p 2?/ - PlyP2x ) ■ 

This can be rewritten as 

(1.12) g e = Ll (\Vp\ 2 + ^\Vr\ 2 ) 

(L 2 + L 3 ) r x 2 r y 2 

+ g ((Piz + y+Jfy) +(P2*-Pii/ + y) ) 

if L 2 + L 3 > 0, 



(1.13) (/ c = (L 1 + L 2 + L3)(|Vp| 2 + ^|Vr| 2 ) 

- g ((y- P1X-P2J,) + (P2* -Ply ~ y) +|Vr| ) 

if > L 2 + L 3 . 

The following structural conditions are assumed for g b (p,t) = <#,(|p| 2 , r): 
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( i)g b € C°°([0, oo) x R),g b > and <? b (^, f ) = 0, 



(1.14V 



ii) For some mi, 7712,7713 > 

IS6,p(|p|V)||p| + Kt(|p|V)| < ™i(|p| 3 + \r\ 3 ) + m 2 , 
m 3 {\p\ 4 + \r\ 4 ) - 1 < g b (\p\ 2 ,r), 



iii) For some S, 777,4 > 

m4({\p\ 2 ~ + \r - l\ 2 ) < g b (\p\ 2 ,r) 
for ||p| _ M| + | r - !| <tf. 



Since f b (Q) = g b (\p\ 2 ,r) = g b (p,x) under the change of variables (1.9), these 
are additional assumptions on f b . From (1.2), (1.12), and (1.13) we see that 
g e is a positive definite quadratic. Thus G £ is strongly elliptic. It follows 
that minimizers for G e in A$ exist and that the Euler-Lagrange equation 
is a semi-linear elliptic system for which minimizers are classical solutions 
(C°°(n) n C 2 (Tt) in our case). (See Theorem 2.2.) 

In general the bulk energy well for ^(|p| 2 ,r) corresponds to {(p, r) : g b (\p\ 2 ,r) = 
0}. From our assumptions on f b , the bulk energy well for f b is U 
{•s( e 3 <g> e3 — |/)}. By the change of variables Q — > (p, r), A' s corresponds 
to T s := {(p, r): |p| 2 = \,r = |}, and s(e3 ® e% — |j) corresponds to 
(p, r) = (0, We note that the structural conditions (1.14) only require 

that {g b = 0} contains T s as in (i), that it is bounded as in (ii), and that g b 
has quadratic growth away from T s as in (iii). 

For the classic example, f b = fP from the liquid crystal literature (with co- 
efficients 0, b, c, and D as described above (see (1.7)), /° minimizes precisely 
on the uniaxial well A s , 



f b {Q) = a(2|p| 2 + |r 2 )-2br(|p| 2 -^) 

+ |(2|p| 2 + ^ 2 ) 2 + c,=:^(|p| 2 ,r) 

for Q e S and Q = Q(p,r), and one can easily show that the structural 
assumptions (1.14) are satisfied for this example of g b . 
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1.2 Main Results 



In this section we state our main results on the structure of minimizers of 
the energy functional F £ (Q) over Aq, using the fact that Q is a minimizer 
of F e in Ao if and only if (p, r) is a minimizer of G t in Aq and Q = Q(p, r). 

Theorem A. Let {(pj,rj)} be a sequence of minimizers for{G £j }, respectively 
over Aq such that Ej 10. For ease of notation we consider pj as a complex- 
valued function by identifying M? and C Then for a subsequence {(pf , rj>)} 
there exists a harmonic function h G C 2 (f2) and k points {a±, . . . , a^} C Q 
such that 

(1.15) (|p / (x)|,r J /(x)) -»• (— ,-) inC loc (n\{ai,...,a k }), and 




in w iw(®\{ a ii ■ • ■ > a k}) n C Joc (n\{ai, . . . , a k }) and in C,™(0\{ai, . . . , a fc }) 
/or a// m > 0, where 0£ = 9e(x) denotes the polar angle of x with respect 
to the center ag. In particular, for each sufficiently small p > 0, if f is 
sufficiently large, setting Q p = Q\ U#=i B p (ae), we have 

(1.16) Pj'(x) = | P/ (x)|e i( V( x )+£?=i m H p 

where hj>(x.) is a function in C 2 (Q p ) so that e*V( x ) has degree zero on dfl, 
and pj/ has degree 1 about each of the k defects {a±, . . . , a^}. 

From the above result and the change of variables between Aq and Ao, we 
obtain: 

Corollary A. Let {Qj} be a sequence of minimizers of {F £j }, repectively 
over Ao such that £j \, 0. Then for a subsequence of minimizers, we have 
Qj' = Q(pji,rf) where {pj',rj/} C Ao satisfies Theorem A, and hence for 
each sufficiently small p > 0, if f is sufficiently large, we have: 

Qj'(x) = Sj' 1 (x)(m J /(x) <g) ny(x)) + Sj> 2 (x)(mj;(x.) <8> mf,(x)) 
- + Sj> 2 (x))I in n p , 
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where 



m/ (x) = (co S (^(/ i/ (x) + ^^(x))),sin(^(/ l/ (x)+^^(x))),0), 

1=1 1=1 

3 3 
Sj'i(x) = |pjv(x)| + -^-/(x), Sj/ 2 (x) = -ry(x) - |pj/(x)|, 



and Qji has degree | about each ag. (See Figure 1.2.) 



In particular, Qf(x.) converges to a uniaxial field Q*(x) in W^(Q\{ai, . . . , cifc})n 
Ci oc (ft\{ai, . . . ,dk}) and in C™ c (Q\{ai, . . . ,ak}) for all m > as j' — > oo, 
where 

Q*(x) = s(m(x) (g> m(x) — - I) in 0\{ai, . . . , a^} when s > 0, 

and 

Q*(x) = s(m _L (x) ® m _L (x) — - I) in f2\{ai, . . . , a^} when s < 0. 

o 

Here 

(1.17) m(x) = (cos(^ (/ l (x) + ^^(x))),sm(^(/ i (x) + ^^(x))),0) 

1=1 i=l 

for all x in 0\{ai, . . . , a^}. Note that rrij/ and m are discontinuous while 
Qji and Q are continuous on Q p . 

The points {a%, . . . , a^} represent the cross sections of the limiting discli- 
nation lines. We prove in this paper that this set of points minimizes a 
reduced energy W(b) defined for b = {b u . . . , b k ) G n k , which was intro- 
duced by Brezis, Bethuel, and Helein in |BB H] in connection with their 
analysis of minimizing sequences {v e } for the Ginzburg-Landau energy 

(1.18) ^ (v) = i^[|Vv| 2 + ^(l-|v| 2 ) 2 ] 

for v G {w G W 1,2 (Q;R 2 ) : w = po/|po| on <9S7}. (The reduced energy W(b) 
is defined by equation (3.28).) More precisely, we have: 

Theorem B. Let {(j>j,Tj)} be a sequence of minimizers for {G Ej } (or equiva- 
lently, let {Qj} be a sequence of minimizers for {F e . } ) for which (a±, . . . , a^) 
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Figure 1: % degree defect in a nematic texture. 



is a limiting configuration of defects as Ej \, as described in Theorem A. 



Furthermore the reduced energy W(b) for the limiting problem minimizes at 
a and we have 



Here 7 is a fixed constant associated to the energy of each defect core. 

The setting we study here gives a good description of two-dimensional ne- 
matic behavior in flat films and thin layers. Investigations from the physics 
literature of nematic textures in flat and curved surfaces (thin shells) can be 
found in [F] , [LP] , [N] , and [VN] . In |FS] Fatkullin and Slastikov propose and 
investigate a model for two-dimensional nematics (assuming that L\ > 0, 
L2 = L3 = 0, and Q = Q(p, r) is in Aq with r(x) = 0) combining Onsager- 



Then 



Fe.iQj) = G £j (pj,rj) - (L 3 - L 2 + \L 3 + L 2 \) 



4 



lim 
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Maier-Saupe and Landau de Gennes theories. This leads them to analyze a 
variational problem closely related to the Ginzburg-Landau energy (1.18). 

Our last result describes how our work in this paper relates to earlier in- 
vestigations of complex Ginzburg-Landau type functionals having multiply- 
connected energy wells. The closest study in this respect is [HKJ by Han and 
Kim in which they analyze the asymptotic behavior for sequences of minimiz- 
ers to the Chern-Simons-Higgs (CSH) and the Maxwell-Chern-Simons-Higgs 
(MCSH) energies used to model aspects of superconductivity. 

For the (CSH) model one seeks (using our notation) minimizers p £ to 

(1.19) C £ (p) = f/ 2 \Vp\ 2 + e- 2 \p\ 2 (l ~ |P| 2 ) 2 ] 

for p G B = {v G W l > 2 (n;R 2 ): v = p on dU}. Here p G C 3 (dft), 
|po| = 1, and deg po = k > with k G N. 

For the (MCSH) model one seeks minimizers (p £j(? , r £i(J ) to 

(1.20) C e ,,(p, r) = J n i\\Vp\ 2 + q- 2 \Vr\ 2 + |p| V + q 2 (e-\\p\ 2 - 1) + r) 2 ] 

for (p,r) G -Bo x ^o' 2 ^)- ^he following two results are from [HK]: 

i) For fixed e > 0, from any sequence of minimizers for (1.20) with q — > oo 
one can find a subsequence {(p £tqe , r £m )} and a minimizer p e to (1.19) for 
which 

Pe, qe Pe and C eM (p eM ,r £iqi ) -> C £ (p £ ) as q t -> oo. 

ii) For fixed q > 0, from any sequence of minimizers for (1.20) with e — > 
there exists a subsequence {(Pe £ ,g, r ££ig )}, a point a q = (of,...,a|) G 
and a function p* as in (1.15) so that p e ^<j —> Pq in the sense of Theorem A 
as E£ — > . 

The functionals (1.10) and (1-20) are quite different. The bulk energy well 
for C £;q is S 1 x {0} |J{(0, e -1 )} and the second component is eventually 
outside of any bounded set as e — > 0. This is in contrast to the bulk energy 
well for G £ which does not vary with e. The analysis in [HK] is based on this 
feature and it cannot be applied to (1.10). Furthermore the bounds in the 
estimates used to prove ii) diverge as q — > oo and so they cannot be used to 
determine lim( lim C e qe (p e qe , r e qe )) = limC e (p e ) or the limiting behavior 

of minimizers of C e , and this was left open. Our analysis however applies to 
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these issues directly. The same arguments we use to prove Theorems A and 
B give the following result: 

Theorem C. Let {p £ } be a sequence of minimizers for (1.19) such thate — > 0. 
Then there exists a subsequence {p £e }, a point a = (a±, . . . , a^) 6 Q k , and 
a function p* as in (1.15) for which p £e — > p* in the sense of Theorem A. 
Moreover W(-) minimizes at a and 

lim [C £e (p £e ) - irk In — ] = W(a) + fc 7 

l— >oo ££ 

for a fixed constant 7. 

Other related work is given in the papers |KSlj . |KS2j . and |SYj in which 
the authors develop asymptotic properties for the (CSH) energy using V- 
convergence techniques. This approach gives less detailed information than 
in our setting. However, it is not restricted to sequences of minimizers as in 
our case, and the authors apply it to more general energies and scalings. 

Our paper is organized as follows. In Section 2 we prove regularity of min- 
imizers and show that minimizers for G £ in Aq correspond to a family of 
equilibria for F £ in A. In Section 3 we prove Theorems A and B, devel- 
oping the qualitative features of minimizers for G £ . Here we expand on 
investigations of minimizers for the Ginzburg-Landau energy E £ (1.18) done 
by Brezis-Bethuel-Helein, Fanghua Lin, and Struwe. (See [BBH], [LI], and 
[St])- The energies E £ and G £ differ in two respects. The elastic term in 
the energy density for E £ is the Dirichlet energy density, whereas for G £ it 
is a coupled quadratic in VQ. Secondly, the energy well for the bulk energy 
density for E £ is S 1 , while the energy well for G £ is a bounded disconnected 
set containing the ring T s as one of its components. The set T s plays the 
same role as the energy well for E £ . For e small we prove that minimizers 
take their values near T s except for an exceptional set contained in a neigh- 
borhood of k defects (vortices). In order to argue as has been done for E £ 
we must first show that this exceptional set has small measure. The results 
in Section 3 are proved assuming the a priori estimate 

(1-21) e~ 2 [ g b (\p £ \ 2 ,r £ )<M 

Jn 

for some constant M < 00, for the family of equilibria {(p £ ,r £ ): < e < 
£1} that are considered. In Section 4 we prove, using a Pohozaev identity, 
that (1.21) is always satisfied if is a disk and < e < 1. We then 
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use this result to establish (1.21) for the case in which Q is a C 3 bounded 
simply connected domain and {(p £ ,r e )} are minimizers, where E\ depends 
on s, L\, L2, L3, Q, k, and the constants in (1.14), and M depends on these 
terms and in addition on ||po||vK 1 > 2 (ao)- Om* approach for this part is similar 
to one used by del Pino and Felmer in |dPF| in which they established the 
analogue of (1.21) for the simpler energy (1.18). 



2 The Landau - de Gennes Energy 



By definition of f e , we have 

fe(Q) 



+ 



^|VO| 2 + ^±^|divO|» 
L 3 

~TT {Qij,kQik,j Qij,jQik,k) 



where div Q is the column vector whose ith entry is the divergence of the 
ith row of Q, Qijj- The last term in f e is a null-Lagrangian; its integral 
over f2 is constant on 

M = {Q e W 1,2 (fi;R 3x3 ): Q = Qo on dtt} 

and its first variation at any element of M is zero. Set 



f'e(Q) 



L 



2 (L2 + L 3 



/|VQ| 2 + 2 



-|divQ| 2 . 



We say that F £ and F' £ = J n [f' e + e 2 ft] are equivalent since their first 
variations on M. agree, 8vF e (Q) = 5vF' e {Q) where 

5 V F(Q) = DF(Q)[V] := 8 t F(Q + tV) at t = 
for all Q G M and V € Wq' 2 (Q; M 3x3 ). 

For Q £ A we write 



<9(x) 



z 2 (x) ^(x) 2:5 (x) 

_z 4 (x) z 5 (x) -zi(x) - z 3 (x) 



and for Q G ^4o we additionally have 2:4 (x) = 25 (x) = 0. The Euler- 
Lagrange equations for F' E derived by variations in _4(„4o) consist of the five 
(three) equations 8 Zl F' e = for £ = 1, . . . , 5 (£ = 1, 2, 3). 
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We can now show that an equilibrium with respect to variations in Aq is 
also an equilibrium with respect to variations in A. We have: 



Lemma 2.1. Let Q G Aq solve 5 Ze F £ (Q) = for t = 1,2, 3, then S Z4 F e (Q) = 
5 Z5 F e (Q) = as well. 

Proof. Since f b = A(det Q, \Q\ 2 ) it is easy to see that d Z4 f b (Q) = d ZB f b (Q) = 
for Q G S - It follows directly that 5 ZA F' £ {Q) = 8 Z6 F' e {Q) = for any 
Q(x) € A- □ 

Theorem 2.2. For each e > minimizers for F e {Q) in Aq exist and are of 
class C°°{n)r]C 2 {U). 

Proof. Recall that by (1.2), L\ > and L\ + L2 + L3 > 0. We consider two 
cases. 

i) L2 + L3 > 0. From the discussion above we can work with the energy F' e 
instead of F e . Its energy density is the sum of nonnegative terms and f' e is 
a positive definite quadratic in Vz, z = (zi, Z2, 23). The first variation of 
F' £ in „4o results in a semilinear elliptic system of three equations in three 
unknowns. From standard elliptic theory (see [G]) minimizers for F' E in ^4o 
exist, they are weak solutions to the resulting elliptic system and they are 
classical (C°°(fi) f] C 2 {U)). 

ii) > L2 + L3. Let curl Q denote the matrix- valued function whose ith 
row is the curl of the ith row of Q. Then |VQ| 2 — |div Q\ 2 — |curl Q\ 2 = 
(Qij,kQik,j — Qij,jQik,k) is a null Lagrangian. As a result, if we set 

/e (Q) = g ' ^' 2^ ' ^' 

then / e — /" is a null Lagrangian, /" is a positive definite quadratic in Vz, 
and we can argue as in the previous case. □ 

Setting pi = (z± — z%)/2,p2 = Z2, and r = z\ + Z3 then Q £ Ao is given 
in terms of (p, r) 6 ^4o by (1.9). The minimum problem for in is 
recast as the minimum problem for G e as defined in (1.10) in Aq, where g e 
as expressed in (1.12) and (1.13) directly corresponds to f' e and f" in cases 
i) and ii) above, respectively. Moreover we have 

! \ f fn\ (-^3-^2 + 1^3 + ^21)/^. ^ r> r> \ 
9e(P, r ) — JeW) ~. Wij,kWik,j — Hij,jHik,k)- 
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Corollary 2.3. If (p, r) G A) and Q = Q(p, r) i/ien 

G e (p, r) = F e (Q) + (L 3 - L 2 + |L 3 + L 2 |) — . 



Proof. It suffices to evaluate J n (Qij t kQik,j — Qij,jQik,k)- As this is a null- 
Lagrangian we are free to choose (p, r) G ^4o, an d we set r = |. It follows 
from (1.9) that 

/ (Qij,kQik,j — Qij,jQik,k) = -4 / (pixP2y - PlyP2x) 

Jn Jn 

= -4fc|B.(0)| = -s 2 7rk. 

□ 

Corollary 2.4- Minimizers (p £ ,r e ) /or G e in Aq exist, they are of class 
C°°(Q) f] C 2 (Q), and they correspond to minimizers for F £ in Aq by the 
relation (1.9). 



3 The Asymptotic Problem 

By Theorem 2.2, equations (1.12)-(1.13), and our assumptions on Qo, it 
follows that minimizers (p e ,r e ) for G £ in Aq are classical solutions to the 
boundary value problem 



(3.1) 



A(p,r) 
^2(p,r) 
£s(p,r) 



= -2LiApi - (L 2 + L 3 )[A Pl + l(r xx -r yy )]=- 2 ^g biP 
= -2L!Ap 2 - (L 2 + L 3 )[A P2 + r XJ/ ] = -^^ p 
= - § Li Ar - \pi xx -piyy + 2p 2xy + \ Ar] = -4^ 

in Q, 



(3.2) and r = -, P = Po on d£l, 

with | po| = ^ on 90 and deg po = k > 0. 

Choose a finite covering of the C 3 manifold by coordinate neighborhoods 
with uniformly bounded C 3 structure, and a constant Eq in (0, 1) (depending 
only on Q. and W) such that for all io £ fi, -B2e ( x o) is contained in a set 
in £/. Throughout this section we assume (1.21) holds for all minimizers 
z e = (p e ,r e ) for G e in ^4q for all < e < e\ < eq, where e\ depends only 
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on s,Li,L2,Ls,Q,k, and the constants in (1.14), and M depends on these 
terms and in addition on llpollw^fdfi)- This will be proved in Section 4. 

We begin this section by proving several a priori estimates, namely Lemma 
3.1 to Lemma 3.6, for solutions to (3.1) and (3.2) that satisfy (1.21) for the 
above M and < e < s±. These and Proposition 3.7 to Corollary 3.13 will 
be applied to minimizers of G £ to prove Theorems A and B at the end of 
this section. 

In this section, unless otherwise stated, we denote by C and Cj positive 
constants depending at most on po, s, L\, £3, £2, and the constants in 
(1.14). Additional dependence, e.g. on M, will be denoted by C(M). 

Lemma 3.1. Let z e = (p e ,r e ) satisfy (1.21), (3.1), and (3.2) for < e < 
E\. Then |z e | and |eVz e | are uniformly bounded in by a constant C(M) 
independent of e for all < e < £\. 

Proof. Let x G H and let e G (0,£i). Set 

z(y) = z £ (ey + x) for y G Cl = {y : ey + x G 0}. 

Then in fi, z satisfies the system obtained by setting e = 1 in (3.1). Let 
B r = B r (0) nfl. From (1.21) and the growth estimate (1.14) on we have 

II^IIl 4 (_Bi) — C(M) for < £ < £1. 

Write (3.1) as £z = £~ 2 f(z), where f(z) = [-2pig bjP , -2p 2 g b)P , -g^xY and C 
is the second order elliptic operator with constant coefficients. From (1.21) 
and the L 4 estimate, we have 

f |f (z) ■ z| < C(M) for <£<£i. 

JB! 

In addition we have 11^1^^^ < C£ for < £ < £\ and £ < 3, where eg 
depends only on and po- 

We use ip 2 (z — tp) as a test function in (3.1) where ip is a cutoff function 
vanishing near |y| = 1, such that ip = 1 on B3/4, and tp is a smooth func- 
tion equal to z on d£l. The above inequalities and elliptic estimates give 
||z|| 12 .£} 34 ^ C(M). This implies that f(z) G L 2 (i? 3 / 4 ) and we see that 

||z|| 2 2-s — C(M). Elliptic estimates imply that z G W^' 2 (Bg/i§) and by 
differentiating the equation we obtain ||z|| 3 2-b 9/16 — C(M). It follows that 
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||z|Lw-g < C(M) uniformly for < e < e±. The assertions then follow 

^ (£>l/2) 

by scaling back to z £ (x). □ 

Set O^: = {(p,r): ||p| - M| + | r - || < M }. Note that O = T s . Below 
% n {E) denotes the n-dimensional Hausdorff measure of E. 

Lemma 3.2. Let z £ satisfy (1.21), (3.1), and (3.2). Set B(s,fi) = {x £ 
f2: z £ (x) £" 0/i} 7 Pi(x,y) = x, and P2(x,y) = y. Let < fi < 5 where 5 is 
given in (1.1 4). Then 

V}{Pi{B{e,n)) < C(n,M)e andU l {P 2 {B{e^)) < C(n,M)e 

for all < £ < e\ . 

Proof. Note that z £ (x) £ Oq for each x £ <9H Let (x',y') £ B(e,[i), and 
set l x i = {(x',y): y £ M}. Since this line intersects d£l there must exist 
(x',y") £ t x i so that z(x',y") £ dO^/2- It follows from Lemma 3.1 that 
there is a Ci(M) > so that 

z(a/,y) £ 3m /4\O m /4 for |y - y"\ < C x e. 

From (1.14) then we see that there exists C<i{p) > so that C2£ < f e t 9b{\Pe\ 2 , f £ )d'H l {y). 
Thus 

C 2 en\Pi{B{e,n)) < [ ^(|p £ | 2 ,r e ) < e 2 M. 
Jn 

The estimate for P2(B(e, fi)) follows in the same manner. □ 
Since B(e,fi) C Pi(B(e,fi)) x P 2 (B(e, n)) for jjl > we have the following. 

Corollary 3.3. Let z £ satisfy (1.21), (3.1), and (3.2). For any [i £ (0,5) if 
< e < e x then H 2 (B(e, /x)) < C(n, M)e 2 . 

This estimate leads to a statement for all x £ J). We use the fact that 
(p £ ,r £ ) is bounded together with Corollary 3.3 for x £ B(e,fi), and the 
growth estimate (1.14) for x £ Q\B(e,fi) to get 

Corollary 3.4. Let z £ satisfy (1.21), (3.1), and (3.2). If < e < E\ then 
(3-3) a' 2 j[((r e (x) - |) 2 + (| Pe (x)| 2 - j) 2 ) < C(M). 
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Lemma 3.5. Let z £ satisfy (1.21), (3.1), and (3.2). If < e < e± then 

[ |Vr e | 2 < C(M). 
Jn 

Proof. We first record an energy estimate for linear elliptic systems applied 
to (3.1) and (3.2), 

l|Pe||i,2;n + Il r £ll2,2;n < c l ( e_4 ( II Pe9b,p + 1 1 ^b,tr 1 1 1;0 ) + IlPO |||,2;dn) 

where c\ depends on L\, L2, L3 and Q. Since gi minimizes on Oq we have 
(3.4) | 56iP (|p £ (x)| 2 ,r £ (x))| 2 + |^ r (|p £ (x)| 2 ,r £ (x))| 2 

<C((|p £ (x)| 2 -^) 2 + (r £ (x)-|) 2 ). 

Thus using (3.3) we find 

lks|ll,2^ < C(^- 2 + l). 

It then follows from this inequality and (3.3) that 

/ |Vr £ | 2 = - / (r £ - ^)Ar £ < £ _1 ||r £ - ^||2;Q£||r £ || 2 ,2 ; n 
Jn Jn <3 o 

< C(M). 

□ 

Lemma 3.6. There is a constant £2 & (0, £1] depending only on and 
k = deg po, and a constant C{M) independent of e so that if (p £ ,r £ ) is a 
minimizer for G £ in Aq and < e < £2 then 

r s 2 1 

/ |Vp £ | 2 < — 2tt kin - + C(M). 
Jn 4 £ 

Proof. We first construct a comparison function for the energy in (1.18). 
Choose a set of distinct points {61, . . . , &&} C f2, depending only on Q, and 
k such that 

min{|6 n — be\, dist(6 ra , dfl); 1 < n, I < k, n 7^ £} = e 
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is maximal. Define 

where ((t) G C 2 {R) such that C(*) = for t < \, ((t) = 1 for 1 < t, and 
j £ (-) is harmonic in 17 such that w e = ^ on 90 for e < e. Then one has 
^e(w £ ) < 7rfcln(^) + co for < e < e where E £ is given in (1.18) and cq 

depends only on Q, and po- We next set (w',r') = (^w £ , |) G Aq and use 
this as our comparison function for G £ . Set £2 = niin {e, ei}. Then for 
£ G (0,£2]j using (1.11) and (1.14) we find that 

G £ (W,r>)<(L 1 + ^±^) [ |Vw'| 2 
z Jn 

+|Li + L 3 | / (w[ w' 2 -w[ w' 2 ) + C\. 
Jn 

The second integral on the right depends only on w'\gQ. Thus we get 

G £ (w',r') < (L 1 + L2 + Ls ) ^ 2vrfcln - + d. 

2 4 e 

Next we use 

( g e (Vp £ ,Vr £ ) < G £ (p £ ,r £ ) < G £ (w',r'). 
Jn 

From (1.11) and suppressing the subscript e we see 
IT , (L 2 + L 3 ) , f |T7 |2 , (L 2 + L 3 ) /, 

(£1 H ^ ) / ' Vp l H 9 / ~ r:rf,2s/ r?/P2a; ^ 

f s 2 1 

+ |L 2 + -L 3 | / (pi x P2y - PlyP2x) < + L 2 + L 3 ) — Trk hl(-) + Cl- 

4 £ 

Again the third integral is a constant depending on po- The lemma will 
follow once we show that we can bound the second integral appropriately. 
To do this we multiply the third equation in (3.1) by (r — |) and integrate 
over Q. We get using Lemma 3.5 that for < e < e 2 : 

{PIxTx ~ PlyTy + P2xTy + P2yT x ) \ 

< £ - 2 Jjg b)K \.\ r - S -\+C 2 (M) 

< e~ 2 [ (\g b ^ + \ r - S -\ 2 ) + C 2 {M). 
Jn 6 
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Finally using (3.3) and (3.4) we see that the last integral is bounded by a 
constant C(M) independent of e for < e < £2- □ 

We are in a position to apply Lin's Structure Proposition, see [L3j . Signif- 
icant parts of the proposition were also proved by Sandier [5] and Jerrard 
Pj. Define 



•M v ) = / ie( v )> where 
Jn 



Proposition 3. 7. For fixed s / and a constant K suppose that 

p £ G {v G W 1>2 {£1) R 2 ) : v = p on dQ} such that 

po € W^ 2 (dn), |po| = |, de 5 (p ) = ft > 0, 

s 2 1 
J £ (p £ ) < vr — fcln - + K 

where < e < n. Fix < a® < min(|, 2 (fc+i) )- There are positive constants 
rjo G (0,??) and po depending on K, Q,po, «o so i/iai if e < t]q then for 
each p £ there are points {a\ , . . . , a|} /or which 

min{|a^ — a||, dist{a £ n , dQ); I < n, £ < k,n ^ £} > po 

and constants a m {e), «o < a m < 2ao /or 1 < m < fc so £/ia£ |p £ | > ^ on 
<9-B m and deg (p e \dB m ) = 1 where B m : = B £am (a £ m ). 

Moreover 

(3.5) / j £ (p £ )< ^ -A- lnj + d 

/or some constant c\ = c\(K, Q, po)- 

Furthermore for any sequence {p £f } wit/i £f | £/iere exists a subsequence 
i £ e(q)}> points {ai, . . . , a^}, and a function h(x) so that 



aZ M -+ a m and p Et(q) -+ p* = {[ M 



. x - a r , 

m=l 
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as q — > oo where the convergence is strongly in L 2 (Q), weakly in 

W loc (^\{ a i> ■ • • > a k}), and \\h\\ w i,2(n) < C2 
for some constant C2 = C2(K, J7, po). 

We take into account (1.21), Lemma 3.5, Lemma 3.6 and apply the Propo- 
sition to a sequence of minimizers. 

Lemma 3.8. Let {(p £ ,r £ )} be a sequence of minimizers for {G £ } in Aq such 
that e X 0. T/ien for a subsequence {(p £e ,r £e )} we have p £e — > p* as in 
Proposition 3. 7 and 

r £i ^ S -m W l > 2 (n). 
The next two lemmas strengthen the notion of convergence using the fact 

k 

that we are working with a sequence of minimizers. Set Vt p = Q\ \J B p (aj). 

Lemma 3.9. Let {(p£,r^)} be a sequence of minimizers for {G £l } in Aq 
converging to (p*, §), in L 2 {Q) where p*(x) = W \[ k j=l j^J e^ x ). T/ien 
/or eac/i < p < ^, 

(p*,r<) -»• (p*, J) in VF 1,2 (n p ) and lime^ 2 / o b (|p^| 2 , r £ ) = 0. 

Moreover Ah = in $7. 

Proof. By Lemma 3.8 and an argument by contradiction we have (pe, rg) 
(p*, |) in M /1 ' 2 (r2 p ) for each p > as above and (p£,r^) is a local minimizer 
for 

/ [« e (Vp,Vr)+e7 2 « 6 (|p| 2 ,r)]. 

To prove strong convergence it is enough to show that for each x £ f2\{ai, . . . , a^} 
there exists a neighborhood of x , on which (p^, rg) — > (p*, |) in 
W 1,2 (U^ n fi). We first consider the case x <9f2 and take d = d(x) > 
_ _ 

such that = -B 2 ^(x) C f2\{ai, . . . , a&}. Then Yl #j( x ) + M x ) is single 

j'=i 

valued here and we write p* = ^ e M x ) on S 2 ^. From Lemma 3.8 there 
exists Cb(M) < oo independent of £ so that 

||p/||l,2;B 23 + lk<||l,2;B aS < C (M). 
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Thus for any subsequence {(p^. ,r^)} of {(p^r^)} (possibly after passing to 
a further subsequence that we do not relabel) d can be chosen, d < d < 2d 
so that 

(3.6) / [|9 r p,.| 2 + \d T r e / + ej*((\p e / - y) 2 + (rt. - S -) 2 )} < d(M) 

JdB d «J 

where 5 r denotes the tangential derivative. Thus (|p^. |,r^ ) — > f ) uni- 
formly on dB d and (p^,r^.) (p*, §) in W l ' 2 (dB d ). Since deg p*\dB d = 
it follows that deg p^-|as d = for j sufficiently large and we can write 

P£.(x) = |p^.(x)|e 1 ^^ for x 6 S-B^. We define tD^.(x) and £>(x) on as 
the harmonic extensions of uj£ j \QB d and w|aB d respectively. It follows that 

(3.7) Q t . ^ w in W l ' 2 {dB d ) and -> w in W 1,2 (S«i). 

The first limit follows from [HKLJ and the second follows from elliptic reg- 
ularity theory. We next construct comparison functions 

(Pij > hj ) ■ = ( I Plj I e l£J ^ , hj ) on £ d 
such that (ptpf^) = (p^,r^.) on &B d . 
This is done by setting 

(Ip^I,^.) = (y, |) on Bd_ e ^, 

and for each define (|p£.|,f^.)(|x|,0) to be linear for d — eg j < |x| < d. 

Then based on (3.6) and (3.7) it follows that (p^ , r t .) -)• (p, f) = e**, f ) 
m W 1 ' 2 Moreover 

/ 5 e ( Vp , 0) = lim / [ff e ( Vp* • , Vf £ . ) + e7 2 g b ( | p* . | 2 , f/, ) ] . 

JB d JB d 3 

From the minimality of (pi j , r£. ) and the weak lower semicontinuity of J B g e 
we have 

(3.8) / <7e(Vp*,0) < limsup / [g e (Vp £j ,Vr £] ) + ej 2 g b (\p ej \ 2 ,r e] )} 

JB d j^oo JB d 

< ( <? e (Vp,0). 

JB d 

From (1.11) it follows that J B g e (Vp,0) minimizes m the set {p = ^ e if G 
W x ' 2 (B d ): f = u> on <9£ d } if and only if A/ = in B d . Thus p is the 
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unique minimizer and p = p* on B&. From (1.12) and (1.13) we see that 
f B ^g e ( 1 \/p, Vr) is the sum of weakly lower semi-continuous integrals. We 
have shown that the sum is weakly continuous on the sequence {(p^ , r £j)}- 
It follows that each of its terms is weakly continuous on this sequence as 
well. Thus j B \Vp £j \ 2 -> f B \Vp*\ 2 and J B |Vr^| 2 -»• as j — > oo. Thus 
(p£-,r^.) — > (p*, |) in W 1 ' 2 (Bd) and as a result the full sequence (pt,re) — > 
(p*, |) in VF 1,2 (i?^). A further consequence is that 

lim ej 2 / gb(\pe\ 2 ,re) = 0. 

a 

Moreover we have shown that p* = ^ e l (Ej=i j+ /l ( x )) w here A/i = in 
n\{ai, . . . ,afc}. From Proposition 3.7 we have that h E lF 1,2 (f2) and this 
implies that the singularities are removable. 

Lastly if x E <9f2 we take a neighborhood and d E (0, eo) so that there 
exists a smooth diffeomorphism defined on B^ satisfying = x and 

^: B+ = {y + x: y 2 + y 2 < d, y 2 > 0} -> 

onto 

We can then carry out the radial construction of (|p^|,f^) in £Tf, push this 
forward to U^, and then argue as in the previous case. □ 

We next prove that {(|p^|,r^)} converges uniformly to (^, |) outside of a 
neighborhood of {a±, . . . ,0^}. The proof is similar to that in [Llj Theorem 
A. This is possible since the density g e can be expressed as the positive 
definite quadratic (1.12) or (1.13). 

Lemma 3.10. Let (p £e ,r £e ) = (pe,rt) be a convergent sequence of minimizer s 
for {G £e } in Aq as in Lemma 3.9. Given p E (0, ^) and fi E (0,5) there 
exists £q so that 

(pi(x), r^(x)) E O u for all x E O p and I > £q. 

Proof. Assume there exists x^ E Q, p such that 

{p t {^), n {^)) £0 M for I EN. 

We consider two cases, 

i) dist(x £ ,dft) > ef for all I, 

ii) dist(x*,dfi) < e°° for all £, 
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where we fix < ao < mm(|, 2 (k+2) ) f rom Proposition 3.7. 

We treat case i) first. Based on (3.3), Lemma 3.5, and Lemma 3.6 we can 
select ctQ < a£ < 2ao and set B' £ = B £ <* e (x^) so that 



2 6 \2 



•s 



2 



n2 



(3.9) / [(|9 rP ,| 2 + \d T r e \ 2 ) + ^(((p, 

+\rt- S f)]<C 
for a fixed constant Cq(M). Define 

for y € Bi(0) = B x . It follows that 

(3.10) |Vp£| + |VrJ| < Cief - 1 , 
for Ci = Ci(M) and (3.9) becomes 

(3.H) / [(l^| 2 + |a T r;| 2 )+ £ f- 1 )((|p^ 

JdB 1 

+1^ - || 2 )] < Co- 

We can assume that B'^ C ^ p /2 f° r each t. Then using |Llj Lemma 1 for the 
first inequality and (3.5) for the last we find 

|degp^| 9Bl | — (1 - a e ) In — - C 2 < / j' (i-« 4 )(p*) 
4 £t JB 1 

f . / \ S 2 7T k , 1 „ 

7b; 4 fc + 1 ££ 

Since < it follows that deg p^|aBi = if £ is sufficiently large. 
Moreover (p«,r») is a local minimizer for 

/ [9e(Vp,Vr) + ef^- l) g b (\p\ 2 ,r)]. 

We can then construct comparison functions just as in Lemma 3.9, and these 
lead as in the previous proof to 

(3.12) lim ef"^ I g b (\v' t \\ r' t ) = 0. 

€->cjo J Bl 
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On the other hand, using (3.11) it follows that 

(p/(y)>ri(y)) G °^ for a11 y G dB ± and * ^ 

By hypothesis, we have (p^(0), r^(0)) 0^. Thus there exists G i?i with 
(p'g(zi), r^((zf)) G 5C 3ai / 4 . Using (1.14), (3.10), and assuming n < 5 (where 
5 is from (1.14)) we see that there are positive constants Cd,f3, depending 
in addition on \x and M, so that 

ffb(|pK x )lVK x )) > /3 for x € B^d-c,^^). 

Thus we conclude that 

sf^ [ g b (\p' £ \ 2 ,r' e )>C 5 >0 

for a constant C§(n, M) > and all t sufficiently large. This contradicts 
(3.12). 

In case ii) we consider (p e< , r e< ) on i? 3e a o (y £ ) p $7 for G <9$7 with |x^— y^| < 
2e"°. We can then flatten the boundary to construct comparison functions 
as in the previous lemma. □ 

In the next two lemmas we prove that if a sequence of minimizers {(p ££ , r e , )} 
converges in (f2\{di, . . . , a^}) then in fact it is bounded in 
W^ c {p\{ai, . . . , afc}) for all j. Our arguments are based on three features, 
first that {(|PeJ>T"ef)} converges uniformly to (^, |) on K for each K CC 

2 

Cl\{ai, . . . , afc}, second that (^, | ) is a nondegenerate minimum point for 
gb, and third that g e is strongly elliptic. A corresponding result is proved 
for minimizing sequences to the Ginzburg-Landau energy (1.18) in [BBH . 
In that case the Euler-Lagrange equations are diagonal and the authors are 
able to apply estimates for elliptic equations. Here our arguments rely only 
on L 2 estimates for elliptic systems. 

Lemma 3.11. Let {(p £e ,r £e )} be a sequence of minimizers for {G £e } in Aq 
converging in W^(Q\{ai, . . . , a^}) as — > 0. Then for K CC fi\{ai, . . . , a^} 
there exist constants £q and E so that if £ > Iq then 

\\D 2 (p £e ,r £e )\\ 2 . K <E 
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Proof. It suffices to establish the estimate in a neighborhood of each point 
in f2\{ai, . . . , a^}. We first consider the case of x$ 6 f2\{ai, . . . , a^}. Then 
-B2d(xo) C f2\{ai, . . . , cik] for some d(xo) G (0, £o). Fixing xo and 77, < 
rj < we take d and £0 so that 

(3-13) / (\Dp £ f + |/> £ J 2 ) < n 

J_B 2 d(xo) 

and 

(3-14) ||p e J - y I + Ve t ~ || < V on S 2d (x ) 

for all I > £ . 

Let C £ C 2 (-B2d(xo)) be such that £ = 1 on £?d(xo). We suppress the sub- 
scripts and write (p ££ , r £ J = (p, r). Then multiplying (3.1) by —d Xj (( 2 d Xj (p, r 
we get using the strong ellipticity of the system that there exists a constant 
A(Li, L 2 , L3) > for which 

A\\(Dd X] (p,r)\\l B2d +ei 2 [ ( 2 [D 2 g b ](d Xj (|p| 2 , r)) • (fty(IPlV)) 

<C|||^C|^.(p,r)||I. B2d - £ 7 2 / 2 56iP |^p| 2 C 2 . 

Here 2?<?{, = (<9p<ft,, d t gb) and [D 2 ^] is the Hessian of gj,. Using (1.14), (3.14) 
and taking rj sufficiently small we have 

Xf C 2 |^(|p| 2 ,r)| 2 < I C 2 [^]^(|p| 2 ,r)-^(|p| 2 ,r) 
J Bid J Bid 

for some A > 0. 

From equations (3.1), using |p| > ^ on Bid, we get 

W ( 2 (9t P + 9l)=ei 4 [ e\Vg p \ 2 <c! C 2 |^ 2 (p, 0| 2 . 

J Bid J B 2d J B 2d 

Thus we find 

(3.15) ||C^ 2 (p,r)||i iflM + e^KV 9b \\l B2d + ^^(IpI V)!!^ 

< Co [ ( 2 \Dp\ 4 + d 

J B 2d 

< Cj t 2 \D 2 p\ 2 .f \Dp\ 2 + C 3 . 

J Bid J B 2 d 
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The last estimate follows by applying the Sobolev estimate 

(3.16) (/ v 2 ) 1/2 <c! (|£M + M) 

with ip = £|Dp| 2 and c = c(f2). Choosing 77 small in (3.13) the first term on 
the right of (3.15) can be absorbed into the left and the lemma is proved for 
the case of K CC fi\{ai, . . . , a^}. 

Assume next that xq G d£l and d < Eq, so that B 2( i(xo) is contained in a 
coordinate patch in which we can locally flatten dQ near xq. We consider 
the special case where dQ is already locally flat, 

B 2d (x ) n (fi\{ai, . . . , a k }) = 

B td( x o) = {(^i) x i) '■ ( x i ~ x oi) 2 + (x2 ~ x 02 ) 2 < 4:d 2 and x 2 > x 02 }. 

Let C e C c °°(S 2( i(a;o)) such that C = 1 on B d (x ). Let (p,f) G W 2 > 2 {VL) such 
that (p, f) = (po, §) on 3f2. Again suppressing subscripts, we multiply (3.1) 
by d Xl (( 2 d Xl (p — p, r — f)) and integrate by parts. Then for any < < 1 
we get 

(3.17) A\\(Dd Xl (p,r)\\l B+ < Ci|||2?C||Q B1 (p,r)|||^ 
+ eeJ%CVg b \\l B+ +\{f |^ lP | 4 C 2 + C 2 ). 

We next multiply (3.1) by 

-d X2 (C 2 d X2 (p,r)) = -( 2 d 2 X2 (p,r) - 2(d X2 (d X2 (p,r). 
Using the ellipticity of C we get 

(3-18) ^IK 2 C(p,r)||^ + -A 1 (||C 2 DQ B1 (p,r)||^ 

+ \\\D(\\D(p,r)\\\l B+ ) 

< - f £(p,r)-d X2 (( 2 d X2 (p,r))=I 

where Ai = Ai(Li, L 2 , L3). 
From (3.1) we have 

1 = [2pi£r 6) p,2p 2 S&,)3,5&,t]* -d X2 (( 2 d X2 (p,r)). 
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Here we integrate by parts. Since g b minimizes at (|p| 2 ,r) = (s 2 ,|), it 
follows that <7b iP = g bjt = on dfl. Thus the boundary term will vanish and 
we find that 

(3.19) L = -ef [ d X2 [2 Pl g b>p ,2p 2 g btP ,g btX ] t ( 2 d X2 (p,r) 
< 2ef I \g b , P \\d X2 p\ 2 C, 2 . 

J Bt d 

Combining (3.17), (3.18), and (3.19) we see that there exists A2(Li, L 2 , L3) > 
so that 

A 2 (\\C 2 D 2 (p,r)\\l B+ + e 7 4 ||CP<7 b ||^ + )<C 2 |||^CI| J D(p,r)||| 2 . B+ 
+ ee^\\CVg b \\l B+ +h[ \Dp\\ 2 + C 3 ). 

From this point the argument proceeds just as above. In the general case 
one first flattens the boundary and analyzes the system in local coordinates 
in the same manner. □ 

Lemma 3.12. Let {(p £e ,r £e )} be the sequence of minimizers for {G £e } from 
the previous lemma. For each integer j > 2 and set K CC J)\{ai, . . . ,a&} 
there are constants Ej so that 

Up ei ,r ei )\\j,2;K < Ej fore>£ . 



Proof. Choose 77 < l f so that [V 2 g b ] > XL on O v . We suppress the subscript 
££ and assume that i > £0 where £0 is from the previous lemma. We further 
assume that d G (0, £0) is sufficiently small so that i?d( x o) C fi\{oi, . . . , a^} 
and so that (3.14) holds. Assume that there exists a constant E q < 00 so 
that 



(3.20) ||(p, r) ||2 2;Bd+£7 2|| ( | p |2_^ )r _^||2_ i2;Bd 
+ £i 4 \\(9b,p, gb,r)\\a)-2,2;B d < E q 

holds for q = j — 1. We prove this estimate for q = j where Ej—i is replaced 
by a possibly larger constant, Ej and d by d/2. Note that we already have 
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(3.20) for q = 2 from Lemma 3.11. Let d 1 be a derivative of order j — 1 and 
L* 9 be the collection of all partial derivatives of order q. Let ( € C'^ ) {B ( i) be 
such that C = 1 on We use (— l) j,_1 3 7 (C 2 9 7 (p, r)) as a test function 
in (3.1) and find 

(3.21) A||CW(p,r)|||^ < C||pC|^(p,r)||| n 

- e,- 2 / c 2 ^(^ P 2 P ,<7 b , r ) • a 7 (p,r) = i - n 

From (3.20) we have / < C (^-i, eZ). We write 

(3.22) 3^, p 2p, ^ >t ) • ST(p, r) = Sn{g h ^3b,x) ■ (2p ■ d 7 P, 3 7 r) 
+ ^ a a d a g b>p dPp ■ <9 7 p, 

\a\<j-2 

(3.23) 2p • a 7 p = d 7 |p| 2 + J] badPp-dfp, 

l<M<j-2 

and 

(3.24) ^(<7 b , p , ff6it ) = [^ 2 <76]^(|p| 2 ,r) 

+ E c ^ II (5 Q |P| 2 )^- II ^r) m ^ 

J2 a ,^W a + \P\mp)=j-l \a\<j-2 \P\<j-2 

where a a ,b a are constants, Iq = tuq = 0, and c a p(x.) = (c^(x), c^(x)) are 
bounded. Inserting (3.22), (3.23), and (3.24) into the right side of (3.21) we 
have for B d = B d (x ): 

II = ef [ C 2 [^]^(|p| 2 ,r)-^(|p| 2 ,r) 

JB d 

+£ t 2 I C 2 Vc a/3 n(a a |p| 2 )^(n^rr^^(| P | 2 ,r) 

JB d ^ 

+ef [ C 2 ^6 )P (E 6 ^ a P-^P) 

JB d 

= III + IV + V + VI. 
Just as in Lemma 3.11 we have 

Xef [ C 2 |^(| P |V)| 2 <///. 

JB d 
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From Sobolev's theorem the derivatives in IV of order less than j — 2 are 
bounded. It follows then for any 9 > that 

i-2 



\IV\ < Cref [ C 2 (El^dP| 2 ' r )| 2 )l a7 (IP| 2 ' r )l 

JB d t=l 

< Betl C A \D^{\v\\r)t + ^f^ d) 
Then using (3.16) and (3.20) we see 



|2 \|2 , C4(Ej-i,d) 



\iv\ < ec z {E^)ef ! c^'-HQpIV)! 2 + 

JB r , 



<B d 

To estimate \V\ we write d 1 = d x iffl' for some x' and integrate by parts to 
get 

-4 / a2i nJ-2 |2 /i^Y /ci \ / a2i tv? |2 C§{Ej-\^d) 
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\V\ <6ej* [ C 2 \D j - 2 9b , p \ 2 + OCsiE^) f ( 2 \D>p\ 2 + 

To bound \VI\ we first consider the terms with a^O. For these |/3| < j — 1 
and we see we can bound these terms just as was done for V . The term with 

2 

a = can be bounded by ^j-^pr + 6\(D : >~ 1 p\ 4 '. The integral of the first term 

over B d is bounded from (3.20) and the second by 0C 8 (-Ej-i) f B ( 2 \D j p\ 2 + 
C 9 (£7j_i,d). Thus 

< c 10 (^--i)(?( e 7 4 jf c 2 i^- 2 ^i 2 + jT C 2 I^Pl 2 ) + C7ll( ^- 1 ' d) - 

Summing on | | = j — 1 and collecting the estimates for III , . . . , V/ we find 
(3.25) A / C 2 |^(P,0| 2 + A(£7 2 / C 2 \D j ~\\p\ 2 ,r)\ 2 

J B d J B d 

< OC^EjM [ ( 2 \DJ(p,r)\ 2 + ef [ C 2 \D^ l (\p\ 2 ,r)\ 2 

V JB d JB d 
-4 f a2|tv?-2/ n m2^ i CtfjEj-^d) 



From (3.1) we have e E 2 {g btP , g btX ) = -(^ ■ (A, £ 2 ), £ 3 )(p, r). Using this, 
the estimate |p| > ^, and Sobolev's theorem we get 



[ C 2 |^'- 2 (56,p,ff6,r)| 2 < Ci 4 (^--l) / C 2 |^'(P,r)| 2 + C 15 (^-l,d). 
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Inserting this estimate into (3.25) and choosing 6 sufficiently small we obtain 
(3.20) for q = j and d replaced by d/2. □ 

Corollary 3.13. Let {(p £e ,r £e )} be a sequence of minimizers for {G £e } in Aq 
converging to (p*, r*) in W^(0,\{a±, . . . , a k })- Then for each integer m 

(Pe e ,r £e ) ->• (p*,r*) in C Joc (n\{ai, . . .,a k }) 

and in C ; ™ c (0\{ai, . . . , a k }) as i ^ oo. 

Proof of Theorem A. Let {(p e ,r e )} be a sequence of minimizers for {G e } 
in Ao for which (1.21) holds and such that e 1 0. Then by applying 
Lemma 3.8 it follows that there exists a subsequence {(p ££ , r e< )} and points 
{ai, . . . , afe} C 0, so that 

(p...^)-(^n^^ , .2)=(P-,i) 

mW^(n\{a 1 ,...,a k })xW 1 ' 2 (n). 

By Lemma 3.10 for each p G (0,£o)> (|PeJ> r e*) — ► (^>|) uniformly on 
_ _ k 

£l p = Q\ |J Bp(aj), and from Lemma 3.9 

(p e „r e J->(p*,|) in W 1 ' 2 ^). 

Moreover h(x) is harmonic in Q. 

Finally, by applying Corollary 3.13 we see that 

(p £e ,r £e ) -> (p*, |) in C(H p ) and C^((l p ) 

for each integer m. □ 

We need to establish several properties for the following minimum problem 
in order to prove Theorem B. Let (3 £ C, \(3\ = 1 and define 

L(-;P) :=L(-,l;P) = L(e,n;P) 
li /i 

(3.26) = inf I [g e (Vv,Vr) + e- 2 g b (\v\ 2 ,r)} 

(v,r)e2l (9 J Bfi 

\s\ 2 e 
+ (2L 1 +L 2 + L 3 ) U vrln(-) 
4 (i 
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where 

13 \ S \ X S 

%p = {( v ,r) G W 1 ' 2 (B fJi ) : v(x) = ^ — and r(x) = - for \x\ = fj,}. 

A \X\ O 

Lemma 3.14- L(t; f3) is independent of {5 for all (3 G C wii/i = 1. 
Moreover L{r) := L(t; (3) is a nondecreasing function of r for r > swc/i 
i/ia£ 7 := lim £(t) > — oo. 

tJ,0 

Proof. For any T G S'0(2), consider the change of variables by rotation, 
y = Tx for x G B\ and set 



in *i2 0' 

*21 *22 
1 



reads as po(x) = yKi for 



The energy density is frame indifferent and as such satisfies 

/e(V y Q(y)) +r- 2 f b (Q(y)) = / e (V x Q(x)) +r - 2 A(Q(x)) 

where Q(y) = i?Q(T*y)-R*- This translates into a statement of invariance 
for g e and 5-5, 

5e(V y p(y), V y f (y))+T- 2 56 (|p(y)| 2 , f(y)) = <? e (V xP (x), V x r(x))+r- 2 56 (|p(x)| 

where p(y) = T 2 p(T*y) and f(y) = r(T*y)- Let = ft + i/3 2 . Then the 
boundary condition for p(x) as a vector in 
|x| = 1 where 

[ fa fa . 

Given T G 6*0(2) the boundary condition for p(y) becomes po 

(y) = M T 2^ T t y 

for |y| = 1. In particular if we let T = K l we get Po(y) = ^y for |y| = 1. 
Thus the mapping (p, r) G 21/3 — > (p, f) G 21 1 is an isometry such that 
G T (p,r) = G T (p,r). In particular we see that L(r;/3) = L(r; 1) = L(t). 

The monotonicity property of L{t) follows by the same argument for (1.18) 
given in [BBH], Chapter 3. A lower bound m for minimizers for the energy 
(1.18) with 0, = B\ is proved in [BBH], Chapter 5. Let u £ be such a 
minimizer with u e (x) = A on dB\. If (v T , r T ) is a minimizer for (3.26) with 
fi = 1 and e = r it follows that 



EJ— v T ) > ^r(u T ) > -vrln(r) 



m. 
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Thus using (3.3) we have 

- / |VvJ 2 > 7rln(r) — m . 

2 J Bl ' Tl " 4 V ; - 

The existence of a finite lower bound for L(r) follows from this and the 
estimates in the proof of Lemma 3.6. □ 

Proof of Theorem B. The relation between F e and G e is proved in Corollary 
2.3. We establish the asymptotic relation by arguing as in |BBHJ, Chapter 
8. Let 

T = {b= (b,...,b k ) G n k : bi^bj if i ^j} 
and for b € T set 

q 6 (x) = h n ^ hb(x) 

j=l 3 

where h b (x) is harmonic in and is determined (mod 2tt) by the condition 
= Po ° n dQ. From |BBHj . Chapter 8 we have 



(3-27) \ [ k |Vq b | 2 = 

1 Jn\ u B p (b 3 ) 

3=1 

s 2 1 

— (vrA; In - + W(b)) + O(p) as p -»• 



where W(b) is the renormalized energy for (1.18) given in [BBHJ . We express 



this using our notation. Set i?(x) = Y2j=i ^ n l x ~~ bj\ an d t = v 1 - where v is 



the exterior unit normal to dQ. Then 



(3.28) W(b) = -^Y\og\h-b J \ + \f R8 U R 

7T? 1 Jan 



&3 

h h d T R+l [ \Vh h \ 2 . 
an 1 Jq 



Note that using (1.11) we have 

a 

\L*2 + Ls\(qbl, x Qb2,y — Qbl,y<lb2,x) 



(3.29) ff e (Vq b ,0) = (Xi + ^^^ lVqbl 2 
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and that qbi, x <lb2, y - qhi, y qb2,x = since |q b | = V- 

We next construct a comparison function for (1.10). Let b E T. Then for 
< ££ « p and for p sufficiently small (depending on Q and b) we define 

r (q b ,s/3) forxeO\U-=i B^bj), 
(Pe e ,r £e ) = < (v (x - bj), s/3) for p/2 < |x - bj\ < p, 
( ( Vi ((x - bj)),rj(x - bj)) for x € B p/2 (bj). 

Here (v„ ri ) minimizes f [de+^ 2 9b] with boundary conditions ft, $; 

on dB p / 2 (0) and /3j = f\ % 3 -b \ e ihb ^\ The function vo is a minimal har- 

monic map valued in {|v| = |||} such that p £l is continuous. From Lemma 
3.14 we have 

(3.30) f [&e(Vv i ,Vr J -) + e7 2 56(|v i | 2 ,r J -)] 

Jb p/2 (o) 

2 

= (2L 1 + L 2 + L 3 )^ln( o ^)+7 + 0e (l) 

as si — >• 0. Then from (3.27), (3.29), and Lemma 3.14 we get 

G £e (p £t ,,r £e ) < G(p £e ,f £e ) 

s 2 1 

= (2L 1 + L 2 + L 3 )-(TTkln(-) + W(b)) + k 7 

4 Eg 

+ 0(p) + o £ (l). 

Let a G Y be a limiting configuration as in Theorem A. Then from Lemma 
3.9 and (3.27-30) we have 

s 2 1 

G £e (p £e ,r £e ) > (2L 1 + L 2 + L 3 )-(7Tk\n(-) + W( a )) + k 7 

4 Eg 

+ 0(p) + £ (l). 

Just as in [BBH]. choosing Eg = £i(p) « p with p — > we arrive at our 
assertion. It follows from these two inequalities that W minimizes at b = a 
and that the limit for G £e (p £e , r £l ) as £ — > oo is established. □ 



4 The Pohozaev Identity 

In this section we show that (1.21) always holds for minimizers of G £ in Aq if 
fi is simply connected and < e < z\ where E\ depends on s, L\, L 2 , L3, O, k, 
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and the constants in (1.14), and M depends on these terms and HpolliyWgn) 
as well. We first prove (1.21) for solutions to (3.1-2) in the case of a disk 
using the Pohozaev identity. 

Lemma 4-1- Let (p, r) = (p £ ,r e ) be a solution to (3.1-2) where 0, = £Ir = 
.Br(O) and < e < 1. Then there is a constant Co = Co(R, Li, £2, £3, ||po||i,2;9B fl) s) 
so that 

-2 



/ 9b(\p\ 2 ,r)<C . 
Jb r 



Proof. We multiply the system (3.1) by — V(pi,]?2, r)x and integrate over 
B R . We find 

(4.1) 0= f [(2L 1 + L 2 + L3)(Ap.Vp-x) + (^ + ^t^)ArVr-x 

Jb r 2 4 

- £- 2 V<7-x] 

+ t / P^Vpa • x + 2p 2:r3/ Vr • x] 
1 JB R 

+ ^ 2 ± L3 ^ /" [( raw - r yy )Vpi ■ x + (pi^ - Pi yy )Vr ■ x] 

* J B r 

= . I+ (L2 + L3) Jf+ (^ + ^) fJJ . 

2 2 
We can calculate I as in [BBH] . Chapter 3, 

(4.2) / = R(L 1 + ^±^1) [ (|p,| 2 -| Pr | 2 ) 

1 J 8Br 

+ R{ ^ + i I ^±M ) f (M2 -| rr | 2 ) + 2£ - 2 / 9b . 

Here p T and r T are tangential derivatives. Note that r T = and p T = po r 
on dBp;. To calculate II we write 

/ r xy SJp 2 • x = / (r xy xp 2x + r xy yp 2y ) 

J Br J Br 

= - (xr x p 2xy + yr y P2 xy ) 
Jbr 

+ p / xy(p 2x r x +P2 y ry). 

K JdB n 



36 



Using this and the fact that r T = on OBr we get 

2 



II 



R 



JdB R 



xyp2 v r v . 



To calculate III we change variables, x' = (x — y)j a/2, y' = (x + y)/ y/2. 
Then 

III = 2 I (r x , y ,Vpi ■ x + p lx , , Vr ■ x) = — / x'y'p lv r v . 
JBr K JdB R 

Writing (x,y) = (RcosO, Rsin9) then it follows that (x',y') = (Rcos(9 + 
f),i?sin(0+f)). Thus 11 + III = R / flB r„(cos20,siii20) -p,,. Finally we 
see that 



(4.3) 



L2 + £3 



(// + III) 



< R 



\L 2 + L 3 \ 



Thus using (4.1), (4.2) and (4.3) with (1.2) we get 



R[L 1 + 



(L 2 + L 3 ) 



iPOrl 



> R(L 1 + 



dB R 

L2 + I>3 \L2 + L3 



\Pu\ 



dB R 



'3Li (L 2 + L 3 ) |L 2 + L 3 | 

1 fl, -r + — s r~ 



dB f 



J Br J Be 



□ 

Lemma 4-2. Let Q, be a C 3 bounded simply connected domain in M?. There 
is a constant < E\ < £0 such that if (p, r) = (p £ ,r £ ) is a minimizer for G e 
in Aq and < e < e±, then 



.-2 



/ g b (\p\ 2 ,r)<M. 
Jn 



Here e\ depends on s, L±, L 2 , £3, k, and the constants in (1.14) a nd M 
depends on these terms and IIpoIIvf 1 ' 2 ^)- 

Proof. Set R = 2(diam(f2)) and assume that £ 0,. We construct an 
extension of p. Let p G W 1,2 (Br(0)\Q) valued in {|p| = ^} and such that 
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p is a minimal harmonic map satisfying p = po on <9S1 and p(x) = 
on dB R (0). Note that \\p\\i,2;B R {o)\n < C||Po||i,2 ; an- Set 



(p, r) for x£H, 
(p,f) forx€£fl\fi. 



Let G £ = f [g e + <#,], and let (p,f) be a minimizer for G £ such that 
Jb r 2e 

(p, f) = (p, |) on <9-Br. We can apply Lemma 4.1 (with e replaced by \ple) 
and the results from Section 3 to G £ and (p, f) for the case of = -Br. 
In particular from the proof of Theorem B there are constants C\ and 
< r\\ < 1, depending on s, Li, L 2 , L 3 , SI, /c, and the constants in (1.14) so 
that 

s 2 1 

(2Li + L 2 + L 3 )j 7rHn--Gi<G £ (p,f)<G £ (p',r') 
for all < e < r/i. Note that 

G £ (pV) = ^[ 5e (Vp,Vr) + ^ «7 b (|p| 2 ,r)] 
+ / 5e(Vp,0) 

= G £ (p,r)-^ J^(\v\\r) + C 2 
where G2 depends only on ||po||i,2;9n and the constants in (1.14). Thus 
(4.4) (2Li + L 2 + L 3 )j7rfclni + ^2 J 9b{\p\\r) < G £ (p, r) + d + G 2 . 

Next we consider the comparison map (w', r') constructed in Lemma 3.6 
defined for e < e = 772- Since (p, r) is a minimizer for G £ we get 

s 2 1 

G £ (p, r) < G £ (w', r') < (2Li + L 2 + L 3 ) jvrfc In - + G 3 

for all e < e = r\\ where G3 depends only on po, ft, Li, L 2 , £3, and the 
constants in (1.14). It follows from this and (4.4) that 

e- 2 / g b (\ P \ 2 ,r) < 2(Ci + G 2 + G 3 ) = : M. 
Jn 

for all < e < ei = min{?7i, 772, £o}- 

□ 
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